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ABSTRACT 

We present a method of constraining the properties of the y-ray emitting region in flat spec¬ 
trum radio quasars (FSRQs) in the one-zone proton synchrotron model, where the y-rays are 
produced by synchrotron radiation of relativistic protons. We show that for low enough val¬ 
ues of the Doppler factor <5, the emission from the electromagnetic (EM) cascade which is 
initiated by the internal absorption of high-energy photons from photohadronic interactions 
may exceed the observed ~GeV flux. We use that effect to derive an absolute lower limit of e); 
first, an analytical one, in the asymptotic limit where the external radiation from the broad line 
region (BLR) is negligible, and then a numerical one in the more general case that includes 
BLR radiation. As its energy density in the emission region depends on 5 and the region’s dis¬ 
tance from the galactic center, we use the EM cascade to determine a minimum distance for 
each value of 5. We complement the EM cascade constraint with one derived from variability 
arguments and apply our method to the FSRQ 3C 273. We find that 6 > 18-20 for B < 30 G 
and ~day timescale variability; the emission region is located outside the BLR, namely at 
r > IORblr ~ 3 pc; the model requires at pc-scale distances stronger magnetic fields than 
those inferred from core shift observations; while the jet power exceeds by at least one order 
of magnitude the accretion power. In short, our results disfavour the proton synchrotron model 
for the FSRQ3C 273. 

Key words: astroparticle physics - radiation mechanisms: non-thermal - galaxies: active: 
individual: 3C 273 


1 INTRODUCTION 


Blazars are a class of Active Galactic Nuclei (AGN), whose broad¬ 
band photon spectrum is dominated by non-thermal emission. This 
is believed to be produced within a relativ istic jet oriented at a small 
angle with respect to the line of sight iBlandford & Reesl 1 19781 : 
llJrrv & Padovanill 19951 ). The spectral energy distribution (SED) of 
blazars is comprised of two broad non-thermal components: a low- 
energy one, that extends from the radio up to the UV or X-ray fre¬ 
quency range, an d a high-energy one that covers the X-ray and y- 
ray energy bands (Padovan i & Giommilll995l : lFossati et alj ri998). 

It is commonly believed that the low-energy blazar emis¬ 
sion is the result of electron synchrotron radiation, with the 
peak frequency reflecting the maximum e nergy at which elec¬ 
trons can be accelerated (e.g. IGiommi et alj 120120 . However, the 
origin of their high-energy emission has not been yet settled. 
Among the proposed mechanisms for y-ray production in bl azars 
are: s ynchrotron self-Compton radiation fe.g. iM araschi et alJ 
Il992l : iBloom & MarscheJ 1 19961 : iMastichiadis & Kirk 1997t) . ex- 
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temal Compton scattering (e.g. 

Dermer et al. 1992|; Isikora et al. 

1994; Ghisellini & Madaul 1996 

), proton synchrotron radiation 

(Aharonian 2000 

; Miicke & Protheroe 2001), and photohadronic 

interactions (e.g. 

Mannheim & Biermann 1993; Atovan & Dermer 

20oT 1; Petropoulou et al. 20151). For Flat Spectrum Radio Quasars 


(FSRQs) in particular, which are characterized by large values 
of the so-called “Compton dominance”, i.e. large ratios of the 
peak high-energy luminosity to the low-energy one, the SSC sce¬ 
nario is disfavoured, w hile the EC and pro t on synchrotron scenar- 
ios remain viable (e.g. Isikora et alj 120091 : IChatteriee et ai] ]2013l: 
iBottcher et al.ll2013h . Besides the radiative process responsible for 
the blazar high-energy emission, the distance of the emission re¬ 
gion from the super-massive black hole that lies in the galac- 
tic center (sub-pc vs. p c scale), remains a ma t ter of debate (e.g. 


Blazeiowski et al. 20001; Tavecchio et al. 201C 

; Poutanen & Stern 

20lB; Marscher et al. 20121; Stern & Poutanen 

2014). In the case 


of FSRQs, the presence of external photon fields may be used to 
constrain the location of the y-ray emission region, at least within 
the leptonic EC scenario (see e.g. lNalewaiko et alJl2014l) . 


In leptonic models the external radiation field has a primary 
role in the formation of the SED, as it provides the seeds for in¬ 
verse Compton scattering in the y-ray regime. In leptohadronic 
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models though, the role of external photons in producing the ob¬ 
served SED is not straightforward. Besides the internally produced 
low-energy synchrotron photons, external photons, e.g. from the 
Broad Line Region (BLR), act as additional targets for photo- 
hadronic interactions with the accelerated protons. In particular, 
if the emission region is located within the BLR, its energy den¬ 
sity as measured in the respective comoving frame will appear 
boosted, thus increasing the efficiency of photopion production. It 
is noteworthy that FSRQs have been suggested as promising sites 


is noteworthy that FMCQs nave been suggested as promts 
of PeV neutrino emission (e.g. lAtovan & Dermerl 2001; 


120131: iMurase et alj I20l4 lANTARES Co llaboratio n et al 


Stecker 

20151 ) . 


Given the recent detection of high-energy astroph ysical neutrinos 
dlceCube Co llaboration! I20I3I : lAartsen et alj l2014h . leptohadronic 
models pose an attractive alternative to leptonic scenarios for blazar 
emission (e.g. lHalzen & Zall997h . 

In this study, we present a method for constraining the Doppler 
factor and the location of the high-energy emission region in FS¬ 
RQs within the proton synchrotron scenario (f or con straints in the 
leptonic scenario of blazar emission, see e.g. iDondi & Ghisellinil 
1 1 99,4 iRani et alj 120131: iNalewaiko etafl l2014l : IZachariasl 120151) . 


Our method is based on the effects of the unavoidable, additional 
emission produced through photohadronic interactions, namely via 
Bethe-Heitler pair production and photopion production. In addi¬ 
tion to the y-ray photons produced by neutral pion (jt°) decay, 
both channels of photohadronic interactions lead to the injection 
of highly relativistic electron/positron pair j], that will contribute to 
the photon spectrum as well. Depending on the parameters that de¬ 
scribe the emission region, such as its size and the magnetic field 
strength, secondary pairs lose energy preferentially through syn¬ 
chrotron or inverse Compton processes, and their emission signa¬ 
tures may appear on the SED (see e.g. lPetropoulou & Mastichiadia 
I2OI5I) . For the case of FSRQs, their emission emerges typically 
at energies much higher than the peak of the high-energy compo¬ 
nent of the SED. and is also_subjected to intrinsic photon-photon 
absorption l!Dernier et ~ai]l2007l) . If the optical depth for intrinsic 
photon-photon absorption (r yy ) is much larger than unity, then an 
electromagnetic (EM) cascade will be initiated, transferring energy 
to the GeV-TeV energy range; it may even dominate over the pro¬ 
ton synchrotron emission (e.g lMannheim et alii 199 ll ) and, in such 
a case, the SED may no longer resemble that of a typical FSRQ. 

An important quantity in our study is the photon/particle com¬ 
pactness, which is a dimensionless measure of the photon/particle 
energy density. It is usually expressed as l oc L/R, where L and 
R are, respectively, the comoving luminosity and size of the emis¬ 
sion region. Roughly speaking, higher compactnesses of the inter¬ 
nally produced synchrotron photons, of the external photons, and 
of the primary injected protons, result in higher photohadroni c pro- 
duction rates and higher optical depths r yy (see also lDermer et alj 
120071 ). This is another manifestation of the so-called “compactness 
problem” in blazars: the photon compactness in the y -ray emission 
region of a blazar cannot become arbitrarily high because of the 
initiated EM cascades that deform its multi-wavelength emissiorQ. 
Since there are different combinations of the Doppler factor and 
of the compactness of primary particles that result in the same ob¬ 
served flux, it follows that the choice of low Doppler factor values 
favours higher production rates of secondary pairs and higher T yy . 


1 From this point on we refer to them commonly as ‘e lectrons’. 

2 This has been also pointed out in jPetropoulou & Mastichiadisli201 2i) for 
the case where no soft photons are initially present in the region (automatic 
photon quenching). 


By reversing the aforementioned argument, it is evident that for a 
given observed y-ray flux, size, magnetic field strength and photon 
compactness, there is a minimum Doppler factor value that can be 
well defined. Taking into account that the total (internal and exter¬ 
nal) photon compactness depends, in turn, on the location of the 
emission region in the blazar jet, we can define, for each Doppler 
factor value, a minimum distance as well. 

This paper is structured as follows. In Sect. [2] we derive an 
analytical expression of the minimum Doppler factor in the limit¬ 
ing case where the internal photon compactness is much larger than 
the external one; this sets the most stringent limit on the Doppler 
factor. In Sect. [3] we present the model and the algorithm for the 
numerical determination of minimum distance of the emission re¬ 
gion. We present the results of our method when applied to the 
FSRQ 3C 273 in Sect. [4] and discuss other possible constraints. We 
continue in Sect. [5] with a discussion of our results and conclude in 
Sect. [6] 


2 ANALYTICAL APPROACH 


The existence of a minimum Doppler factor (<5 m i n ) in the proton 
synchrotron model for blazar emission can be demonstrated with 
analytical arguments. In doing so, we will also derive an analytical 
expression that will reveal its dependence on the quantities describ¬ 
ing the blazar emission region. Our analysis will be focused on the 
minimal case where the internally produced synchrotron photons 
are the only targets for photohadronic interactions; this is also the 
case of minimum photon compactness and may be realizable if the 
emitting region is located much further out of the region of external 
radiation (see Fig. 0. 

We approximate the electron synchrotron differential luminos¬ 
ity as a broken power-law: 


L'(e') = A 1 e'^ 1 H[e’ - e’] + A 2 ^ k H[e' - e'], (1) 


where H[x] = 1 for x > 0 and 0 otherwise. From this point on, we 
use the following convention: primed and unprimed quantities are 
measured in the comoving frame of the emission region and in the 
observer’s frame, respectively. The normalization constants are 


A 2 

■At 


(1-/LX/L-1) , A , , 
02 -Pi £S S 


( 2 ) 

( 3 ) 


where 0 2 > 1, < L L' s = L s /S 4 is the total synchrotron lumi¬ 

nosity, and e' = e s (l + z)/6 is the break energy of the synchrotron 
spectrum, which for the particular choice of spectral indices coin¬ 
cides with the synchrotron peak energy. 

The differential number density of synchrotron photons is 
written as 


ra'(e') = A\e'- r 'H[e[ - e'j + A 2 e'- ri H[e' - e'], 


( 4 ) 


where T, = 0 X + 1 and T 2 = 0 2 + 1 are the low- and high-energy 
photon indices, respectively, and 


- 


AncR 2 


^1,2- 


( 5 ) 


The optical depth for the absorption of y-ray photons with energy 
e[ is 


T yy (e[) ^ 


Rcr 0 (m c c 2 ) 2 


*J2(m e c 2 , 


de'^1 

J 2(m e c 2 ) 2 / ^ 


( 6 ) 
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where we approximated the photon-ph oton absorption cross sec - 
tion aj^] cr yy (xix) ~ crnH\xiX-2V(xix) ( jCotmi & Blandfordll990b . 
Here, xi,x are the photon energies in units of m e c 2 and tr 0 = 
0.652cr T . For e'ej > 2(m e c 2 ) 2 the integral simplifies into 

_ 3<r 0 m,/h) + z) 2 f 

77 1 87 Tfi?< 5 3 (1 + z) e s \ 2 (m e c 2 ) 2 6 2 ) 


where 


(1-#)(#>-1) 

0S 2 -A)(l+y9i)' 


( 8 ) 


The y-rays produced through n° decay are, in principal, very 
energetic and can easily satisfy the threshold criterion for photon- 
photon absorption on the synchrotron photons with energy e'. This 
can be understood as follows. In the proton synchrotron model, 
the high-energy component of the blazar SED is explained as syn¬ 
chrotron radiation of relativistic protons. The maximum proton 
Lorentz factor is related to the peak frequency (v 7 ) of the y-ray 
spectrum as 


Tp.max 


/ 2 /r(l + z)v 7 m p c \ 1/2 

\ q e B6 ) 


(9) 


or, using indicative parameter values, 

Tp.max = 2.4 X 10 s (Bi6i)~ 112 ((1 + z)v %22 ) 1/2 , (10) 

where we introduced the notation q x = q/lO* in cgs units. The 
typical energy of y-ray photons produced by neutral pion decay 
is e^.o^-> 7 — 0.5K pn m v y v cr where k p „ - 0.2 is the mean proton in¬ 
elasticity; in fact, the inelasticity increases from ~ m*/m p (- 0.14) 
close to the thresh old to ~ 0.5 at an energy three time s larger than 
the threshold one llSteckeil 1 9681 : iBetrel man et aT]ll99ol) . Thus, pro¬ 
tons with Lorentz factor y p>max result in the production of very high- 
energy photons: 

<o^ 2y « 2.4 x 10 4 TeV (B^y 112 ((1 + z)v 7 , 22 )‘ /2 . (11) 

For a fiducial synchrotron peak energy e s = 0.1 eV, which corre¬ 
sponds to 


< 


= 0.01 eV 


( 1 +z) 

Si 


e s 

0.1 eV’ 


( 12 ) 


we find that » 2(m e c 2 ) 2 . By substitution of eqs. jilt and 

O in eq. 0 we find the respective optical depth to be 


Tyy(6jfi_^2y) 


4 x 10 3 / O.l eV \ 1/2 L s ,45 

( 1 +z) \ / ^ 16^1 


1+zf 4 

Si I 


Vy,22 

— 


1/4 


(13) 


where we assumed/L = 1/2 and/3 2 = 3/2. Since Ty y (e„ o_, 2y ) 3> 1 is 
typical, the y-ray luminosity from n° decay (L n o _, 2 ) will be totally 
absorbed. We may thus write that Z/*® = (1 -e~ T r>)L n o_, 2y =s 

L jr o^ 2y - The absorbed photon luminosity will be re-distributed at 
lower y-ray energies through the development of an EM cascade. 
This emerges as an additional emission that should be below the 
proton synchrotron component, which in our framework is respon¬ 
sible for the FSRQ high-energy emission. 

We note that photons emitted by secondary, highly relativistic 
electrons from charged (n + ) pion decay or/and Bethe-Heitler pair 
production are also subject to photon-photon absorption, and thus 
they may contribute to the cascade emission. The synchrotron pho¬ 
tons emitted by secondary pairs are less energetic than those from 


3 For simplifying reasons, we neglected the logarithmic dependence . 


7 T° decays. Thus, they are mainly attenuated by photons with ener¬ 
gies e' > e'. To exemplify this, let us consider the most energetic 
pairs produced by pion decays. These are produced roughly with 
y pn si (\IA)K pn y V j iax m v lm e , and the respective synchrotron photon 
energy is written as 


e s,pir — 60 TeV <5j"‘(l + z)v 7i22 . 


(14) 


This lies just above the threshold for yy absorption on photons with 
e'. Cooling of pairs with y pn as well as the production of pairs from 
protons with y p < y p , max results in photon emission at e' <K e' pjr , 
where the threshold condition for absorption on e' is no more satis¬ 
fied. Since n' oc e' -1- * for e' > e', the optical depth for absorption 
of photons with e' <K e' pjr is expected to be much less than 


r yy( e s, P jr) ~ 


10 2 


( 1 +z) 


0.1 eV \ L sA5 
e s ) R^S\ 


(15) 


where we used eqs. 03. ® and €^e' p;r ~ 2(m e c 2 ) 2 . Moreover, the 
synchrotron spectrum from pairs spans many decades in energy and 
the luminosity emitted at e' p;r is, therefore, only a fraction of the 
total injected luminosity in pairs. This is not the case for the y-ray 
spectrum from ji° decays, which is sharply peaked at Simi¬ 

lar arguments apply to the synchrotron emission from Bethe-Heitler 
pairs, which are produced on average with y < y pn . Thus, for the 
purposes of this analytical approach, we can safely ignore the at¬ 
tenuation of photons from Bethe-Heitler and pion decay process, 
and consider only the attenuation of y-ray photons from 7 r° decay. 
In any case, our analytical results will be compared against those 
calculated numerically, after taking into account the additional pho¬ 
ton emission from Bethe-Heitler and charged pion processes (see 
Sect.0. 

Since the proton synchrotron emission alone can explain the 
observed y-ray spectrum, the sum of the cascade and proton syn¬ 
chrotron emission may exceed the observations for high enough 
values of L„o_, 2y . This can be avoided if the following energetic 
constraint is satisfied 


L^ 2y < r, L y . pk , (16) 

where L y pk is the peak luminosity of the high-energy SED compo¬ 
nent, which is typically a good proxy of the total y-ray luminosity. 
Here, r\ < 1 is a dimensionless factor to be defined later by the 
observations (see Sect. 0 The above relation does not take into 
acco unt any spectral information about the dev eloped EM cascade 
(e.g. lMannheiir]| 1 9931 : IPetropoulou et aljfet) 1 3l) . It is based on the 
simplifying assumption that the absorbed luminosity re-emerges at 
the energy where r 77 (e*) ~ 1. Using fiducial parameter values and 
solving eq. 0 for e*, we find 

f * - 18 GeV^j^<5?L- 2 5 R? 6 , (17) 

where we also assumed pS\ = 1/2 and yS 2 = 3/2. We thus expect 
the EM cascade to emerge at energies higher than the peak of the 
high-energy emission, which for FSRQs usuaKyJalls in the 4 MeV- 
40 MeV (10 21 Hz-10 22 Hz) range jFossati et al.lll998h . The con¬ 
straint imposed by relation HI} could be relaxed if we were to in¬ 
clude spectral information for the cascade emission. However, this 
lies out of the scope of the present work. 

The y-ray luminosity from the n° decay may be written as 

L^ 2j ^ L^ 2 y - 2 x l -T pn L p , (18) 

where T pn is the optical depth for photopion interactions (to be de¬ 
fined below) and L p is the total injected proton luminosity. The fac¬ 
tors 2 and 1/2 account for the production of two photons that share 
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the energy of the parent neutral pion. In what follows, we will use 
eqs. d and 4H for deriving and justifying the existence of a 
minimum Doppler factor. 

Before we calculate the optical depth for photopion {pn) in¬ 
teractions, it is useful to determine the threshold photon energy for 
such interactions with protons having Lorentz factor y pmax . Using 
eq. {TOjl we find 

4 “ 0-58 eV (BxSi) 112 ((1 + z)v r>22 )~' /2 , (19) 

which for typical parameter values, is e' h » ^ (see also eq. On. 
Lower energy protons will interact with synchrotron photons of en¬ 
ergy e' > e/ h , whose number density decreases as oc g'-t-fe; we w ip 
not consider these interactions in the following. 

The optical depth for pn interactions is defined as T pn (y p ) = 
t cr / t, m (y„), where f cr = R/c and t pj ) is the pn energy loss rate given 
by Isteckeilfl968l) 


tliji f deecr p „(e)K pn (e) f de'^-^p-, (20) 

2 Tp Jf th Jf/ 2 yp e'- 

where e th = 145 MeV. We assume K p „{e) ~ 0.2 and approximate 
the cross section as cr pn = d- p „H[e- e th ], where & nn = 1.5 x 10 ~ 4 ctt 
( for a more realisti c description of a pn , see iMiicke et al.l l200d : 
iBeringer et alj|20l2h . For the photon spectrum defined by eq. 0. 
the second integral is written as 

de’A 2 (!^- pi . (21) 

For the highest energy protons, i.e. with Lorentz factor y p , max , only 
the second integral is non-zero, and the respective optical depth is 
written as 


r ^ = r + r 

e/2y p € J e/2y p «Jma? 


, . ^ 3d~pjrKpn L s 2 , (3 \ ) /2y p,max + Z) \ 

T^(y p , m ax) - 2ncR ) 


( 22 ) 


where 


g(j32,j8 1 ) = 


(i-Am-i) 


(23) 


p 2 (pi -/?i)(2 + Pi) 

Using eq. dlOt. B\ = 1/2, p 2 = 3/2 and fiducial values for the other 
parameters, the optical depth is written as 


^(Tp.max) — 6 X 10 


L s ,45V: 


,3/4 

y,22 


^i 6 Bf 4 df /4 ( 1+ ^) 1/4 


(—f- 

V0.1 eV/ 


(24) 


Using the approximation L r pk =s L r , the definition of the proton 
compactness C p , which is a dimensionless measure of the proton 
luminosity, given by 


_ <r T L p 
p 4nRm p c 3 8 4 ’ 

and eqs. 1 1 61 . i ll 8b and ( 1221) we find that 


(25) 


6g(p 2 ,Pi)^ A- fed + z)t-- 1 f< p m p c 2 6 < 1. (26) 
o- T r/L y \ e t h / 

At this point, we make use of our working hypothesis, namely 
that the y-ray emission is explained by proton synchrotron ra¬ 
diation. Using standard expressions for the synchrotron luminos- 
ity emitted by a power-law proton distribution (e.g. eq. (6.36) in 
iRvbicki & Lightmaiil i ll 9861) ) we may express { p as 


C p o- T LyB-^ +1),2 6^ p+5y2 

4nR 2 f pV ?- p)l2 


(27) 


where / p = (p - 2)/(p - 1) with p t 1 being the power-law index 


of the injected proton distribution and 

(3 - p){p + 1) (2n) <p ~ ly2 m ( p p+,)/2 c (p+3)/2 

P ” 2ql 5+p)l2 3 p ' 2 r 1 (p)r 2 (p) 



( 28 ) 

(29) 

(30) 


where T(/) = dxx!~ l e~\ The c onstant C p * is a generalization 
of C 2 given by eq. (9) in IPetropoulou & Mastichiadisl J2OI2I) for 
2 < p < 3. Substitution of eqs. (Hi and ( 1271 ) into eq. l !26b results in 


8 > <5min. 


(31) 


with 


A (p+3+3&)/2 _ „ . -l„-2 J32-1 v «J 2 +p-3)/2 r -032+P+1)/2 

°min “ « t s V O 


(32) 


In the above. 


C 1 y(/l 2 ,/l] )(TpnKpj r ttl p C Cp i _ 

Jy Jp \ ^th 


fl h/2 / 2 \* (2nm p c\ f>212 


<le 


(33) 


and / z s 1 + z- It is noteworthy that <5 min does not depend on the 
y-ray luminosity, while it has a weak dependence on most of other 
model parameters: 

x R-l/w jj-2/w p2/wy(fi2 + P ~3)/ w ( 34 ) 


where w = p + 3 + 3/L>. Instead, the strongest dependence comes 
through the magnetic field strength and the spectral index of the 
synchrotron spectrum above its peak: 


• the minimum Doppler factor decreases for stronger magnetic 
fields. Higher values of B require lower proton luminosity to ex¬ 
plain a given observed y-ray luminosity, as eq. {27j demonstrates. 
This subsequently reduces the luminosity produced through pho- 
tohadronic interactions and, thus, the amount of energy being ab¬ 
sorbed and reprocessed (see eq. GU); 

• the minimum Doppler factor decreases as the electron syn¬ 

chrotron spectrum above its peak becomes steeper. This can be eas¬ 
ily understood, since the respective number density of synchrotron 
photons scales as n'(e') oc and decreases for higher p 2 . We 

remind that very high energy y-rays produced via photohadronic 
interactions are mostly absorbed by these synchrotron photons. 


The dependence of <5 min on the magnetic field strength is ex¬ 
emplified in Fig. [0 where <5 m j n is shown as a function of B for two 
values of the emission region radius: R = 3.6 x 10 16 cm (solid line) 
and 3.6 x 10 15 cm (dashed line). Other parameters used for the plot 
are: v s = e s /h = 3.2 x 10 13 Hz, L s = 6.3 x 10 45 erg/s (this cor¬ 
responds to 10 13 Jy Hz for the 3C 273 distance D L = 755 Mpc), 
Vy = 10 22 Hz, p = 0.2, p = 2.3, pi = 0.7 and p 2 = 1.35. Our choice 
of the parameter values is motivated by the SED fitting of 3C 273 
(see Sect.0. 

Figure [7] demonstrates that the Doppler factor of the high- 
energy emission region in FSRQs lies above 15-20, unless the mag¬ 
netic field strength is high, i.e. > 30 G and/or the synchrotron pho¬ 
ton spectrum is steep, e.g. p 2 > 2. It is important to note that we 
have arrived at this conclusion without considering any additional 
constraints imposed by e.g. the observed high-energy variability. 
In other words, even if the observed variability is not faster than 
~hour timescale, the cascade emission initiated by photohadronic 
interactions limits the Doppler factor to values larger than 15-20. 

The minimum Doppler factor shown in Fig. |T| can be consid¬ 
ered as an absolute lower limit, since it was derived using only the 
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B (G) 

Figure 1. Minimum Doppler factor ci nnn as a function of the magnetic field 
strength for R = 3.6 X 10 16 cm (solid line) and R = 3.6 X 10 15 cm (dashed 
line). The region below the curves leads to modification of the proton syn¬ 
chrotron spectrum due to the enhanced cascade emission, and is therefore, 
forbidden. Other parameters used for the plot are: v s = e s /h = 3.2X10 1 ' Hz, 
L s = 6.3 X 10 45 erg/s, v r = 10 22 Hz, T] = 0.2, p = 2.3, fix = 0.7 and 
jS 2 = 1.35. 


internally produced radiation. If we were to include an extra low- 
energy photon component in our calculations, such as the BLR pho¬ 
ton field, the Doppler factor would have to be larger than dmin (see 
Sect.[4]l. 


3 NUMERICAL APPROACH 

In what follows we will expand upon the idea presented in the pre¬ 
vious section by including in our calculations the emission from the 
BLR. Our working framework is analogous to that adopted in the 
previous section but with two main differences: 

(i) the use o f the numerical code described in 
iDimitrakoudis et af] J 2 OI 2 I ) allows us to make no assumptions 
about the photohadronic emission and the initiated EM cascade. 
The steady-state proton, electron and photon distributions are 
self-consistently calculated by numerically solving the system of 
coupled integrodifferential equations that describes their evolution 
in the energy- and time-phase space. 

(ii) the minimum Doppler factor for a particular FSRQ will be 
derived by fitting its multi-wavelength contemporaneous observa¬ 
tions with the proton synchrotron model. 


3.1 Model 

We model the blazar emission region as a spherical homogeneous 
blob of radius R that contains a tangled magnetic field of strength 
B. The region moves with a bulk Lorentz factor T at a small angle 
0 obs with respect to the observer. The respective Doppler factor 8 is 
defined as 8 = T _1 (1 -ft cos (LbsU 1 ■ We assume that relativistic elec¬ 
trons and protons with power-law distributions are being injected 
into the source at a constant rate, while they may physically escape 
at an energy-independent timescale that is set equal to the crossing 
time t a = R/c of the emission region. Electrons lose energy through 
synchrotron radiation and inverse Compton scattering on the exter¬ 
nal photons (EC) as well as on the internally produced synchrotron 


photons (SSC). Protons lose energy by emitting synchrotron ra¬ 
diation and through the photohadronic channel s of Bethe- Heitle r 
pair p roduction and photopion production (e.g.lPimitrakoud is et alJ 
120121) . The loss processes will lead to the injection of secondary 
electrons and photons which are, respectively, subjected to syn¬ 
chrotron/inverse Compton s cattering and photon-photo n absorp¬ 
tion. We refer the reader to iDimitrakoudis et al] l20l4) for a de¬ 
tailed description of the physical processes. 

As already mentioned in the introduction, our working hy¬ 
pothesis is that the low-and high-energy components of the blazar 
SED are the result of primar y electron and proton sync hrotron ra¬ 
diation, respectively (see e.g. iMucke & ProtheroefeOOll ). The emis¬ 
sion produced through photohadronic interactions appears at even 
higher energies than the high-energy component, and is subjected 
to photon-photon absorption. 

We assume that the spectrum of the external ultraviolet (UV) 
radiation arises from an optically thick accretion disk, and is then 
scattered by the BLR clouds. For simplicity, we approximate the 
accretion disk emission with a black-body spectrun|j that peaks at 
th e observed energy. For t he BLR we adopt the geometry presented 
in lNalewaiko et all i j 20141 14 . and illustrated in Fig.0 The inner ra¬ 
dius of the BLR is defined as R blr , whic h is related to the accretion 
disk luminosity L a d iGhisellini & Tavecchidl2008i) as 

Rblr ~ 10 17 cm L*^ 45 . (35) 

Other sources of external photons could be the reprocessed line 
emission or infrared radiation from a dusty torus. In what follows, 
we will not include in our calculations the radiation from the torus, 
since its luminosity and size are less well-defined than those for the 
BLR. We will also ne glect the direct irradiation from the accretion 
disk dPermer & SchI ickeisei] 120021) . T his is a safe assumption as 
long as the emissi on region lies at (e.g. lGhisellin i & Madau] [l996l : 
ISikora et alJl2009ll 


r > 0.8 


D nl 
(TBLR 


1/3 


—■ 0.01 pc 


1/3 L l/3 m 1/3 
BLR.-l ad,45 BH.9’ 


(36) 


where M BH is the black hol e mass, which for 3C 273 is Mrh.q = 
A?bh/(10 9 Mq) - 0.9 - 2.4 dPeterson et alj|2004 : IPaltani & Ttirlerl 
2005), and £ BLR is a dimensionless factor that incorporates all the 
details about the geometry and the irradiation of the BLR from the 
accre tion disk. A representative value is £ B lr ~ 0.1 d Si kora et alJ 
2009|), while values < 0.01 are considered to be very low (e.g. 
Nalewajko etal]l2014l) . 

The emission region, which is depicted as a yellow blob in 
Fig. El is located at a distance r in the blazar jet. In this study, 
we treat r as a free parameter, i.e. r S R blr , with the aim of im¬ 
posing a minimum value on the ratio r/R B lr. As shown in Fig. [2] 
the radius of the emission region is kept constant, with only the 
constraint of being smaller than the transverse size of the jet at a 
distance r, i.e. R < r/T, where we also used 8T < 1 (see 

iNalewaiko et alj|2014l and discussion, therein). We will return to 
this issue in Sect. [5] where we discuss how our results would be 
altered, if we allowed R oc r v . In principle, R is related to the 
observed variability timescale as t var > (1 + z)R/c8. As r var may 


4 An optically thick, geometrically thin disk, i.e. a Shakura-Sunyaev disk, 
is better described by a multi-temperature black body, whose flux scales as 
F(e) oc c 1 /4 exp(—e/eo). However, the details of the accretion disk spectrum 
do not affect our analysis. 

5 Some obs ervations may sug gest, however, a planar geometry for certain 
FSRQs (e.g. lStem & Poutanen|[2014h . 
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gas and dust broad line 

torus Ublr ^ re 9 ion ( BLR ) 



emission region 

(inside and outside BLR region) 



Figure 2. Sketch of the geometry of the emission region and its surrounding 
sources of radiation, for an FSRQ. 


take a wide range of values, even for the same source, depend- 
ing on the observing period and ener gy band (for 3C 273 see e.g. 
iKataoka et al.ll2002| - |Soldi et al.ll2008l) . we choose to use R instead 
of f var as the free parameter. We discuss the variability constraints 
later in Sect. 14.31 

Following Isikora et a 0 d2009h and iNalewaiko et~aD {20141) . 
we write the energy density of the BLR region as measured in the 
rest frame of the emission region as 


0A^ LR r 2 L, d 

3ncR BLR 


m. 


(37) 


where the function A is defined as 


A(x) 


1 

1 + x 4 ’ 


(38) 


and x = v/Rblr- Similarly, the BLR photon energy as mea¬ 
sured in the_comoving frame of the emission region is given by 
jNalewaiko et al.li2014t) 

e BLR - TeBLR-j—• (39) 

An important quantity in our analysis is the so-called photon com¬ 
pactness which is defined as C y = «(,o" T i?/m e c 2 , where u' y is the co¬ 
moving energy density of an arbitrary photon field. Using eq. 03 
we may write the BLR photon compactness as 


(blr - 4r 2 d(x), 


(40) 


where 


4 = 


OA^bLRO - T^4id 


3itm c c 3 R 


2 

BLR 


(41) 


Thus, the total photon compactness in the emission region, which 
is relevant to the calculations of photohadronic emission, is given 
by 4,t = 4lr + dyn, where { syn is the compactness of internally 
produced synchrotron photons. We note that we do not take into 
account the anisotropy of the BLR photon field as seen in the emis- 
sion region in the calculations of photopion p roduction (see also, 
lAtovan & DermeJl200ll : iTavecchio et al.ll2014l) . 


6 Similarly, we have used in Sect. [2] the term proton compactness (see 
eq. {27}). 


3.2 Method 

The algorithm we follow in our numerical approach is described 
below. 


3.2.1 No external radiation 

We start by assuming that 4lr is negligible with respect to the 
compactness of the internally produced photons, namely 4 lr 'K 
4yn- This can be seen as the case of minimal compactness. For a 
given pair of B and R we 

(i) choose a high value for the Doppler factor 6, e.g. 50; 

(ii) choose values for the rest of the parameters, e.g. f p ,c and y Piinax , 
that lead to a reasonable fit of the SED. This is defined by the curve 
that passes within the error bars of most of the observational points, 
with particular emphasis on the highest energy ones, that are more 
directly affected by secondary particles from photohadronic pro¬ 
cesses within this model. We are not interested in the absolute best 
fit, as would be determined by a^ 2 test, but rather a good enough 
one, as determined visually; 

(iii) if the derived photon spectrum describes the SED reasonably 
well, we return to step (i) and choose a smaller value of 6; if the 
photon spectrum does not fit the SED for the adopted Doppler fac¬ 
tor because of enhanced photohadronic emission, we stop and de¬ 
fine the current value of the Doppler factor as <5 min . 


3.2.2 Internal and external radiation 


As a second step, we include the BLR emission into the calculation 
of the broad-band photon spectrum. For each value of <5 ~ T, we 
obtain the multi-wavelength spectra for different values of 4 lr- 
Each pair of (6 , 4lr) translates into a pair of (6 , x) through eq. <40 1 . 
Thus, numerical runs for fixed S and 4lr imply different locations 
of the emission region in the jet, i.e. 



(42) 


This corresponds also to different comoving photon energies e RLR 
(see eq. <39t ). It is important to note that the inclusion of the BLR 
radiation does not affect the values of 4, P and y pmax that we de¬ 
rived previously (Sect. 13.2. It . since the SED is fitted by the syn¬ 
chrotron radiation of primary electrons and protons. This shows the 
secondary role of the external radiation in the proton synchrotron 
model, in contrast to the EC leptonic models. 

We then determine that value of 4 lr above which the cas¬ 
cade emission modifies the proton synchrotron radiation spectrum 
at a few GeV in a way that the total emission exceeds the ob¬ 
servations. As can be evidenced by eq. {42}. the maximum value 
of the BLR compactness translates into a lower limit of r/R blr. 
This parametrization of the problem allows, therefore, for solutions 
within or outside the BLR and is a generalization of the approach 
presented in Sect. l3.2.H We note that the cascade emission depends 
on both 4 lr and e RLR , which affect the photohadronic production 
rates in a direct and indirect way, respectively. The value of e RLR 
affects the energy thresholds for photohadronic interactions. Thus, 
for a given Doppler factor, a choice of a higher value of I BLR does 
not necessarily mean that the luminosity of the cascade emission 
will be higher. 
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4 RESULTS: APPLICATION TO 3C 273 

We apply our method to the well-known FSRQ 3C 273 at redshift 
Z = 0.158. The optical-UV spectrum of 3C 273 shows a prominent 
excess of emission, which is mainly interpreted as a contribution 
of the accretion disk emission (see Ulrich, 1981; Soldi et al., 2008, 
and references therein). The detection of lines in the optical-UV 
spectrum of 3C 273, e g. Ly-or, CIV, OVI, CIII, NIII, and SVI (e.g. 
Paltan i~& Tiirlerl [2003L and references therein) is connected with 
the BLR. The accretion disk and BLR luminosities of 3C 273 are 
well-defined, i.e. L ad = 1.3 x 10 47 erg/s jVasudev an & Fabianl2009t) 
and L B lr = 9.1 x 10 45 erg/s l l Peterson et alJI 20041) . The knowledge 
of both L ad and Lblr reduces the number of free parameters enter¬ 
ing in the model (see also lBottcher et alJl2013l) . Other parameters 
describing the BLR are Rblr ~ 1.1 x 10 18 cm (from eq. 03), 
£blr = 8 eV and £ B lr = 0.1. As can be evidenced by the ASI 
Science Data Centre (ASDC0, there is a huge amount of archival 
and non-simultaneous observations for 3C 273. As simultaneous 
multi-wavelength observations are important for our analysis, we 
use the dataset by lAbdo et ah 1 ( flolch - We emphasize, though, that 
our method can be easily applied to different broad-band simulta¬ 
neous data, since it is based on a generic idea. 

4.1 No external radiation or r » R B lr 

Figure [3] demonstrates the effect of the cascade emission on the 
high-energy part of the spectrum as the Doppler factor of the emis¬ 
sion region progressively decreases, in the minimal scenario where 
only internal radiation is a target for photohadronic interactions. 
The three panels (from top to bottom) correspond to different mag¬ 
netic field strengths, namely B=30 G, 15 G and 7.5 G. The ra¬ 
dius of the emission region is assumed to be R = 3.6 x 10 16 cm 
in all runs. For a fiducial value of 6 = 15, our choic e results in 
t var ~ 1 day, which is typical for 3C 273 ( Icourvoisier et all 19881 ). 
Other parameters used and kept fixed in the numerical runs are: 
ye,min = L6 X 10 , yp,min L ye,max = 5 X 10 , p e = 2.7 and 
p = 2.3; both distributions of primary particles were modelled as 
n, oc y-Pig-T'/ri.max j = e, p. The parameters that had to be adjusted in 
order to model the SED for the different Doppler factor values are 
listed in TableQ] In all cases, the plateau-like emission above a few 
GeV (> 10 24 Hz) is the result of the EM cascade initiated by VHE 
y-rays produced in photohadronic interactions. Spectra shown with 
red thick lines correspond to the minimum value of the Doppler 
factor that can explain the simultaneous SED. For 6 < (5 min , the 
photon spectra above > 10 24 Hz exceed the observations. 

At this point, it is interesting to compare the numerically de¬ 
rived and ci min listed in Table [I] with the respective values pre¬ 
dicted by our analysis in Sect.[2] This is exemplified in Fig.[4]where 
the top and bottom panels show the comparison for C p and <5 m i n , 
respectively. In both panels, the values from the numerical analy¬ 
sis are shown with symbols, while the curves are calculated using 
eq. (|27]l for R = 3.6 x 10 16 cm, v 7 = 10 22 Hz, L 7 = 6.3x 10 46 erg/s 
and p = 2.3 (top panel) and eq. 1321 for v s = e s /h = 3.2 x 10 13 Hz, 
L s = 6.3 x 10 45 erg/s , v 7 = 10 22 Hz, 77 = 0.2, p = 2.3, ft = 0.7, and 
ft = 1.35 (bottom panel). 

In both panels, our analytical curves are in good agreement 
with the numerical values determined through the SED modelling, 
with some deviation becoming systematically larger for B = 7.5 G 
(top panel) and R = 3.6 x 10 15 cm (bottom panel). However, it is 



log v (Hz) 



log v (Hz) 



Figure 3. Proton synchrotron model fits to the simultaneous multi¬ 
wavelength data of 3C 273 by lAbdo et al]|2oi()l (black symbols) for B=30 
G, 15 G and 7.5 G (top to bottom), in the limiting case of £blr Lyn or, 
equivalently r 3> Sblr- In each panel, the SED is modelled using different 
values of the Doppler factor. Spectra shown with red lines correspond to 
the minimum value of the Doppler factor that can explain the simultaneous 
SED. The respective spectrum when proton synchrotron radiation is omit¬ 
ted, which reveals the underlying spectrum of the EM cascade, is shown 
with a red dashed line. The accretion disk spectrum is overplotted with an 
orange dashed line. Grey symbols are non-simultaneous and archival obser¬ 
vations taken from ASDC. 
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Table 1 . Parameter values used for modelling the multi-wavelength emis¬ 
sion of 3C 273, as illustrated in Fig. [3] Not all cases listed in this Table 
are depicted in Fig. [3] but we include them for completeness reasons. Other 
parameters used are kept fixed (see text). 



<5 

4 ( in l°g) 

(p (in log) 

B=30 G 
y P ,max = 10 8 


28 

-5.2 

-3.5 


26 

-5.1 

-3.4 


24 

-4.9 

-3.2 


22 

-4.8 

-3.1 


20 

-4.6 

-3.0 

(^min) 

18 

-4.4 

-2.8 


16 

-4.2 

-2.6 

B=15 G 

7p,max = 1.2 X 10 8 


32 

-5.5 

-3.3 


30 

-5.3 

-3.1 


28 

-5.2 

-3.0 


26 

-5.1 

-2.9 

(^min) 

24 

-4.9 

-2.7 


22 

-4.8 

-2.6 


20 

-4.6 

-2.4 

B=7.5 G 

y p ,max = 1.6 X 10 s 


40 

-6.0 

-3.0 


38 

-5.9 

-2.9 

(Anin) 

36 

-5.7 

-2.7 


34 

-5.6 

-2.6 


32 

-5.5 

-2.5 


remarkable how well the analytical curves follow the trend found 
numerically, and especially for £ m j n , since we made several ap¬ 
proximations in order to derive an analytical expression (eq. [32). 
A quantitative difference between the curves is something to be ex¬ 
pected, since the numerical analysis: (i) takes into account the emis¬ 
sion from secondary pairs (from Bethe - Heitler and n + decays) in 
the formation of the EM cascade, (ii) makes no assumptions about 
the photohadronic production rates, and (iii) takes into account the 
spectral shape of the EM cascade. 

For completeness reasons, we repeated the modelling proce¬ 
dure for a smaller emission region having R = 3.6 x 10 15 cm. The 
bottom panel of Fig. E] shows that unreasonably high values of the 
Doppler factor (<5 > 40) are required, in this case, to avoid the ef¬ 
fects of the EM cascade, while only very strong magnetic fields 
(B » 50 G), can bring the Doppler factor to lower values. 

4.2 Internal and external radiation 

Following the method described in Sect. 13.2.21 we derived the min¬ 
imum value of the ratio r/R B L r for the three cases considered pre¬ 
viously. Our results are summarized in Table [2] and Fig. [5] The 
different curves are obtained for B = 30 G, 15 G and 7.5 G. The 
results obtained from the numerical fitting are shown as symbols, 
while the curves are the result of interpolation. Only the region 
above each curve is allowed, as values below the curve lead to sig¬ 
nificant emission from the EM cascade. Thus, each curve is the 
locus of points corresponding to the minimum distance, r m j„, for 




Figure 4. Comparison between the analytical (curves) and numerical 
(symbols) results. Top panel: t. p as a function of S for B = 7.5 G 
(solid line/squares), 15 G (dotted line/circles) and 30 G (dashed-dotted 
line/triangles). Bottom panel: <5 m ; n as a function of B (bottom panel) for 
R = 3.6 X 10 16 cm (solid line/squares) and R = 3.6 X 10 15 cm (dotted 
line/circles). 

various Doppler factor values. The curves can be also safely ex¬ 
trapolated to higher Doppler factor values, since a power-law de¬ 
pendence is established, namely r m i„/7?BLR 00 <5 1/2 . The abrupt in¬ 
crease of r min /R BLR occurs at 6 = <5 min , as expected. We remind that 
dmm is derived in the limiting case of £ B lr ^ 4yn or, equivalently, 
r » Rblr- 

Figure [5] reveals the following trend: curves move from the 
lower left part to the upper right part of the plot for progressively 
weaker magnetic fields. The horizontal shifting of the curves can 
be easily understood by inspection of eq. <23, which shows that 
d m ; n oc B ( 3+ A)/(p + 3+ -W). For a given S, weaker magnetic fields re¬ 
quire higher values of the proton compactness to explain the ob¬ 
served y-ray luminosity, which in turn enhances the EM cascade 
emission. To avoid an excess in y-rays due to the cascade emis¬ 
sion, the emission region should be located even further out. This 
qualitatively explains the vertical shift of the curves in Fig. [5] 

Figure[5]shows that the emission region of 3C 273 in the pro¬ 
ton synchrotron scenario cannot be located in the BLR region, at 
least for B < 30 G and R = 3.6 x 10 16 cm. Let us discuss how a 
different choice of R and B would affect our conclusion. A choice 
of a smaller radius, e.g. R ~ 10 13 cm, would have a similar effect 
on the curves as that of a decreasing magnetic field (shifting to the 
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Figure 5. Plot of the ratio r/R blr as a function of the Doppler factor for 
three values of the magnetic field marked on the plot and R = 3.6xl0 16 cm. 
Symbols and lines are used for displaying the numerical results and the in¬ 
terpolated values, respectively. Only the region above each curve is allowed, 
as values below the curve lead to significant emission from the EM cascade. 
The abrupt increase of r/Rft lr occurs at 6 = £ m j n defined by eq. {32}. 


Table 2. Maximum /’blr and, equivalently, minimum t//?blr as deter¬ 
mined by modelling of the multi-wavelength emission of 3C 273 for a given 
Doppler factor and magnetic field strength. The method for the determina¬ 
tion or their values is described in Sect. 13421 All other parameters are same 
as in Table Q] 



8 

^BLR 

r/f?BLR 

B=30 G 


28 

10~ 4 

7.5 


26 

10~ 4 

7.3 


24 

10~ 4 

7.0 


22 

10~ 4 

6.7 


20 

10~ 4 

6.4 


18 

1(L 5 

10.8 

B=15 G 


38 

10~ 5 

15.6 


32 

10~ 5 

14.3 


30 

10- 5 

13.9 


28 

10- 5 

13.4 


26 

10- 5 

12.9 


25 

5 x 10~ 6 

15.1 


24 

io- 6 

22.1 

B=7.5 G 


45 

2 x 10^ 6 

25.4 


40 

2 x 10~ 6 

24.0 


38 

2 x lO^ 6 

23.4 


37 

10~ 6 

27.4 


36 

8 x 10~ 7 

28.6 


upper right part of Fig. [5}. Only if the emission region were larger 
and, thus, less compact could it be located within the BLR. How¬ 
ever, as we show in Sect. 14.31 this scenario becomes less plausible 
when the variability of the source is taken into account. In addition, 
the jet power of a larger emission region would significantly exceed 
the accretion power of 3C 273 (see details in Sect. 14.3L 

In principle, the emission region could be located within the 


BLR for sufficiently large magnetic fields, namely B » 30 G, ac¬ 
cording to the trend we find in Fig.[5] The question that arises in this 
case is whether the required B values are plausible or not. Instead 
of performing additional simulations, we can address the question 
with the analytical tools presented in Sect. [2] Inspection of TablesQ] 
and [2] shows that at the minimum distance r min the electron and 
BLR compactnesses are approximately equal, while 8 = d m i n +e, 
with e — 1 - 2. Since the electron compactness is a good proxy for 
the internal synchrotron photon compactness, we can estimate the 
minimum distance by requiring ^BLR(<5 m in) — 4yn(<5mm) or, equiva¬ 
lently 


8 2 . ( 0 

mm u 

1 + 




(43) 


where (5 m j n is given by eq. ( 132b . For R = 3.6 x 10 16 cm and for 
all other parameters same as in Fig. |j] we find that r m m/-RBLR < 1 
for B > 450 G. We can therefore argue that the emission region of 
3C 273 cannot be located within the BLR for plausible parameter 
values. 


4.3 Additional constraints 

The constraints on the Doppler factor and on the distance of the 
emission region from the super-massive black hole were derived 
based only on the photohadronic emission and the initiated EM 
cascade. These can become even more tight when combined with 
information about the variability of the source and the energetics of 
the emission region. 


4.3.1 Variability 


The observed high-energy variability in blazars may range from 
hours up to few days depending on the flaring activity. For example, 
the shortes t variability timescales probed by FERMI-LAT are sev- 
eral hours dTavecchio et alJl2oTol : IFoschini et alJl201 it ISaito et aD 
1201 j) . In order to keep our analysis as generic as possible, we chose 
R and 8 to be independent parameters. Here, we revise our previ¬ 
ous results by including the variability information. If f var is the 
observed variability timescale, 8 and R should satisfy the causality 
condition 


8 > <5 var 


*(l+z) 

Ct V ar 


(44) 


Figure [6] illustrates the revised parameter space 8 - B for R = 
3.6 x 10 16 cm and three indicative values of the observed variabil¬ 
ity timescale marked on the plot. The regions that lie above the 
horizontal dotted and solid lines denote areas where both the EM 
cascade and variability constraints are satisfied. For f var = 1 d, the 
variability offers no additional constraint over 8 min . However, for 
Lar £ 12 h the lower limit 8 rir becomes more constraining than the 
lower limit derived imposed by the EM cascade. 

Similarly, the r/R BLR -S parameter space can be further con¬ 
strained by including the variability constraint. This is illustrated 
in Fig. [7] The vertical lines show d var calculated using eq. ( 1441) for 
R = 3.6 x 10 16 cm and three indicative values of f var (from left 
to right, f var = 24 h, 12 h and 6 h). To the right of those vertical 
lines are regions where the causality condition is satisfied. The EM 
cascade does not exceed the y-ray observations at a few GeV for 
parameters drawn above the curves. Finally, the coloured regions 
denote the parameter space where both requirements are satisfied. 
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Figure 6. Same as Fig.[l]with the addition of the variability constraint 8 ^ 
8 V ar (dotted lines) for R = 3.6x 10 16 cm and three indicative values of t vas 
marked on the plot. The curves of <S m j n are calculated using eq. {32}. Only 
the regions that lie above the horizontal dotted lines, satisfy the variability 
constraint. All other parameters are same as in Fig. |T| 



Figure 7. The revised parameter space t//?blr _ 8 for 3C 273 after tak¬ 
ing into account the variability constraint 8 ^ £ V ar for R = 3.6 x 10 16 cm 
and three values of the magnetic field marked on the plot. Vertical lines de¬ 
note <5 var given by eq. <44V From left to right: t yai = 24 h, 12 h and 6 h. 
The coloured regions above the curves and on the right of the vertical lines 
denote areas where both the variability and EM cascade constraints are sat¬ 
isfied. All other parameters are same as in Fig.Q] 


We find that for ~day timescale variability the location of the emit¬ 
ting region and the Doppler factor are only limited by the EM cas¬ 
cade. A shorter variability timescale, wh ich h as been observed dur¬ 
ing bright y-ray flares of 3C 273 (e.g. Rani e t alj2013h , can impose 
tighter constraints on the minimum r and 8 ; e.g., for f V ar = 6 h and 
B = 7.5 G only the upper right comer of the parameter space is 
allowed. 


4.3.2 Jet power 

In FSRQs the accretion disk luminosity can be estimated using the 
BLR luminosity, and the accretion power (P acc = Me 2 ) can be then 



Figure 8. Jet power, as derived in the proton synchrotron model for 3C 273. 
The values derived by numerically modelling the SED of 3C 273 (see Ta- 
ble|3j are shown as symbols, while the curves are calculated using eq. ED 
for B = 7.5 G (solid line/squares), 15 G (dotted line/circles) and 30 G 
(dashed-dotted line/triangles). The horizontal dashed line marks the accre¬ 
tion luminosity given by eq. {45}. 


calculated as (e.g. lGhisellini et ai]l2014l) 

Pace = 10—, (45) 

< 5,-1 

where e r is the radiative efficiency. In general, the jet power is writ- 
ten as P ]a = fjP a cc, with f) £ 1.5 {Zdziarski & BdttchetfcO 1 si . and 
references th erein). Although the jet po wer can exceed the accre¬ 
tion one (e.g. lTchekhovskov et al.l201 lh due to the efficient extrac - 
tion of energy from a Kerr black hole felandford & Znaiekil 19771) . 
here we consider the more conservative case of 6j = 1. We therefore 
impose the following ‘energetic’ constraint 


p < p 

1 jet ^ 1 acc» 


(46) 


Neglecting the cold proton and radiation energy densities^the 
powe r of a two-sided jet can be written as (e.g. iGhisellini et alj 

s 

Pjet ~ 2nR 2 T 2 c («' + u' p + Kg) (47) 

where «' (i = e, p, B) is the energy density as measured in the rest 
frame of the emission region. Dropping the electron term (see also 
Table0 and using eq. {27} with 8 ~ T, we write the jet power as 

Pje. * Li [ A (Bd)~ (p+1)/2 + (B<5) 2 ], (48) 


where 


A = 


6C p m p c 2 L y 

f P R 2 vt P)l2 


(49) 


The jet power given by eq. {48} for R = 3.6 x 10 16 cm, L r = 
6.3 x 10 46 erg/s, v 7 = 10 22 Hz, p = 2, and three values of the 
magnetic field, i.e. B = 7.5, 15 and 30 G, is plotted as a function of 
8 in Fig. [8] Overplotted with symbols are the values calculated us¬ 
ing eq. {48} for u' p determined by the numerical SED modelling of 
3C 273 (see also Table{3}- Apart from an offset (< 3) between the 
analytical curve and the numerical values for the case of B = 7.5 G 
(see also Fig. 0. the two results are in good agreement. Figure [8] 
shows that in all our SED fits the jet power is dominated by the en¬ 
ergy density of relativistic protons, i.e. Pj a oc (B8T (p+ly2 , while it 
exceeds P icc = 1.3x 10 48 erg/s (dashed line). Even the minimum jet 
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Table 3. Energy densities of relativistic electrons and protons (in the co¬ 
moving frame) as determined by modelling of the multi-wavelength emis¬ 
sion of 3C 273 for a given Doppler factor and magnetic field strength. The 
jet power is listed in the last column of the table. All other parameters are 
same as in Table |T| 



6 

w' a 

u' h 

u p 

P C 
'jet 

B=30 G 

(«g = 35.8 erg/cm 3 ) 


28 

1.7 

58 

1.8 


26 

2.3 

75 

1.8 


24 

3.3 

113 

2.1 


22 

4.4 

150 

2.2 


20 

6.5 

200 

2.3 


18 

11.0 

320 

2.8 


16 

16.8 

502 

3.3 

B=15 G 

(«B = 9 erg/cm 3 ) 


32 

2.8 

100 

2.7 


30 

3.8 

132 

3.1 


28 

5.1 

182 

3.7 


26 

6.8 

245 

4.2 


24 

9.6 

339 

4.9 


22 

13.7 

483 

5.8 


20 

20.5 

753 

7.4 

B=7.5 G 

(u' B = 2.2 erg/cm 3 ) 


40 

2.2 

182 

7.2 


38 

3.0 

244 

8.7 


36 

4.1 

332 

10.6 


34 

5.5 

458 

13.0 


32 

7.5 

621 

15.6 


a Electron energy density in units of 10 -5 erg/cm 3 . 
c Electron energy density in erg/cm 3 . 
b Jet power in units of 10 49 erg/s. 


power (Pjet,min) exceeds the accretion power by approxim ately one 
order of magnitude (see also IZdziarski & Bottchen[2Q15h . Thus, at 
least for the particular choice of R and B, relation < f46t cannot be 
satisfied, and in this regard, it cannot further constrain the parame¬ 
ter space. 

The relatively good agreement between the analytical and nu¬ 
merical results for the jet power allows us to use expression {48} 
for investigating the dependence of ^jet, min on the parameters, and 
to search for those, if any, that can bring the jet power closer to the 
accretion luminosity. The jet power given by eq. {48} is minimized 
for 

I p+ l\ 2/(p+5) 

6 0 B 0 = IA^- (50) 

and its minimum value for a given source (L r , v 7 and p, fixed) de¬ 
pends only the radius R through 

^jecmin = A/f (2 ^ 2 >^ +5 >, (51) 


where 


A = 


6 C p m p c 2 L y 


4/(p+5) 


f yC-ril 2 


-(p+l)/(p+5) + yt/Cp+5) 


(52) 


and 


s = 


ml 

4 


The requirement ^je^min = Pace is satisfied for 


(53) 


Ro = 



(p+5)/(2p-2) 


(54) 


which for 3C 273 and p = 2.3 becomes Rq ~ 10 14 cm. Substitution 
of R 0 into eq. ( 150b results in 8qB 0 ~ 5x 10 4 (in cgs units), namely the 
Doppler factor and the magnetic field should take extreme values. 

Summarizing, we showed explicitly that in the proton syn¬ 
chrotron scenario for the multi-wavelength emission of 3C 273 
there are no reasonable physical parameters that can bring the 
jet power close the accretion power. In all cases, we find 
P\et <: I OPacc. in agreeme nt with the independent analysis by 
IZdziarski & Bdttchei d2015t) . 


5 DISCUSSION 


The derivation of constraints for the Doppler factor and/or the 
location of the emission region in the leptonic framework of 
blazar emission has been the subject of several studies (e.g. 


Dondi & Ghisellini 

1995 

Poutanen & Stern 2010; Tavecchio et al. 

20lT Dotson et al. 

2012 

Cerruti et al. 2013; Dermeretal. 2013; 


Nalewaiko et alJl2014l : IZachariasI 120151) The exploration of such 


constraints within the confines of alternative models is a method 
that may lead to their eventual verification or exclusion, with new 
observations potentially reshaping the available parameter space. 
In this study, we expand this search by adopting the proton syn¬ 
chrotron model for the y-ray blazar emission. 

Under the assumption that the high-energy component of the 
SED in FSRQs is explained in terms of proton synchrotron radi¬ 
ation, we showed that for low enough values of the Doppler fac¬ 
tor (5, the emission from the EM cascade may exceed the observed 
~GeV flux. In fact, the superposition of the EM cascade and proton 
synchrotron components results in a spectral hardening of the total 
emission above a few GeV (see Fig.[3}. From our analysis, it is not 
clear why the proton synchrotron component should dominate in 
the hard X-ray/soft y-ray regime while the EM cascade should be 
suppressed. One could then naturally pose the following question: 
Would it be possible to fit the SED with the cascade component in¬ 
stead of the proton synchrotron component, and would the param¬ 
eters of such a model be reasonable? Roughly speaking, the peak 
of the cascade spectrum can be found by the condition r ry (e*) ~ 1 
(for more details, see lMannheiml dl993i) ). For the parameter values 
used throughout the text, we showed that e* ~ 18 GeV (see eq.ll7t. 
In principle, one could find parameter values that could bring e* 
down to a few MeV, since e* oc L“ 2 /? 2 <5 8 6s. In this scenario, the 
peak luminosity of the cascade should be higher in order to explain 
the observed peak y-ray luminosity, namely r pn L p - L y . 

Unless Tp n > 1, which would correspond to R ~ 10 14 cm, B <; 1 G, 
and 6 < 5 (see eq. {24j ), this scenario would require higher pro¬ 
ton luminosities, and therefore, even more extreme jet powers than 
those listed in Table[3] These estimates, however, are made by con¬ 
sidering only the internal synchrotron photons as targets for both 
photopion interactions and yy absorption. A proper answer to the 
question posed above requires a self-consistent calculation of the 
cascade emission by taking into account the external photon fields 
as targets for both yy and photopion interactions and by focusing 
on a different parameter regime than the one considered here. Such 
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an investigation is interesting on its own, and will be the subject of 
a future study. 

A serious challenge to the proton synchrotron model arises 
from the need for high values of B at distances far from the cen¬ 
tral black hole. lPushkarev et~afl d2012l) have determined, using core 
shift measurements, the magn etic field - distan ce relation to be 
B ~ 0.4 G x (l pc/r) (see also Zdziarski et all2014L for a theoretical 
investigation!. ISavolainen et all d2008h have presented specifically 
for 3C 273 magnetic field measurements of its pc-scale inner jet 
struct ure. A ssuming that the jet’s angle to our line of sight is 10° 
dStawarzll2004h . it would appear that B < 8 G at r - 1.25 pc. Ap¬ 
plying the linear relation between B and r and using R B lr = 0.3 pc, 
we find that for the values of B considered in Fig. [5] and Fig. [7] 
our derived values of r/i?BLR are over a factor of 5 higher than 
those needed to accommodate such strong magnetic fields. This 
would imply that, at least for the case of 3C 273, the magnetic field 
strength required by the proton synchrotron model at the location 
of the emission region is in conflict with the observations. 

In principle, the parameter space used in modelling the SED 
could be further constrained by requiring that the jet power (Pj ct ) 
should not exceed the accretion power (P ac , c ) which, for FSRQs like 
3C 273, can be safely estimated. We showed, however, that in the 
proton synchrotron model for 3C 273 the jet power exceeds that of 
accretion, i.e. ^jet > 10_P acc , for all reasonable parameter values. We 
note, however, that this is as much a problem for leptonic models 
as for hadronic ones dCihisellini et alJl2014l) . 

A possible caveat of our analysis is our assumption of a con¬ 
stant radius for the emission region, as well as of a constant mag¬ 
netic field strength. Assuming that the emission region fills the 
entire cross section of the jet, one could express the size R as a 
function of distance r from the central black hole, i.e. R oc r s 
with s > 0 depending on the specific model of the jet structure 
(e.g. iGhi sellini et~aT] 1 19851 : iModerski et alj 120031 : IPotter & Cotter! 
l2013ll . Similarl y, the magnetic field could be written as B oc P * 
with q < 0 (e.g. lviahakis & K6niglll2004l ; lKomissarov et al]|2007lf . 
Yet, there would remain some arbitrariness regarding the relation 
between the magnetic field in the jet and in the emission region, 
as this would depend on the dissipation mechanism (for a discus¬ 
sion, see Isironi et alJl2015tf , which in turn d epends on the dis¬ 
tance from the supe r -massive black h ole (e.g. ISikora et alj|2005l ; 
iGiannios et alj |2009: Nalewajko 2012). In any case, we can qual¬ 
itatively predict the effects of an increasing size and decreasing 
magnetic field strength on the results presented so far. On the one 
hand, a larger R would loosen the EM cascade constraint, while it 
would push d var to higher values. This could bring the location of 
the emission region closer to or inside the BLR, at the cost of even 
higher jet powers. On the other hand, a decreasing magnetic field 
would make the EM cascade more stringent, i.e. higher <5 mill would 
push the location of the emission region even further away from the 
BLR. The impact on jet power would depend on the Doppler factor 
but it would be marginal, as evidenced by eq. {48}. In short, an in¬ 
creasing radius would have the opposite effects from those of a de¬ 
caying magnetic field. Thus, their combined effect would strongly 
depend on the details of the model, such as q and s. 

The question naturally arises whether the method presented 
here would yield significantly different results when applied to 
other FSRQs, or to a flaring period of 3C 273. Our results regard¬ 
ing the high jet power are not expected to differ, since thi s is a n in¬ 
trinsic feature of the proton synchrotron model (e.g. lBottcher et all 
1 2013l l. It is not straightforward, however, what the answer would 
be regarding the location of the emission region in other luminous 
blazars. The hypothesis that the emission region, in the hadronic 


framework for FSRQs, is located at the ~pc scale jet, can be easily 
tested by applying our method to a sample of densely monitored 
sources, and we plan to do so in a future study. 


6 CONCLUSIONS 

We have demonstrated a method of constraining the properties of 
the y-ray emitting region in FSRQs within the one-zone proton 
synchrotron model. Even though the high-energy component of 
the blazar SED is attributed to synchrotron radiation of relativis¬ 
tic protons, the emission from photohadronic processes cannot be 
avoided. In fact, the EM cascade initiated by the absorption of 
photons produced via photohadronic interactions may exceed the 
observed y-ray flux at ~GeV energies, for small enough Doppler 
factors. For the purposes of our analytical treatment, we focused 
on the photons from neutral pion decay, while we neglected the 
synchrotron emission from Bethe-Heitler pairs and those produced 
through the decay of charged pions. Therefore, the EM cascade re¬ 
sults only from photons from neutral pion decay, which interact 
with the low-energy synchrotron blazar emission. To avoid a sig¬ 
nificant alteration of the SED in that energy range, a constraint is 
set on the luminosity from the EM cascade which is translated to a 
lower limit on the Doppler factor given by eq. {32}. Our analytical 
calculations were performed in the asymptotic limit where the ex¬ 
ternal radiation from the broad line region (BLR) is negligible. In 
this regard, (5 min can be considered as an absolute lower limit. 

We then applied our method to a single FSRQ, namely 3C 273, 
and generalized our analytical method by including the radiation 
from the BLR, as well as the variability timescale (t var ) in our cal¬ 
culations. Because the BLR energy density in the comoving frame 
of the emission region is highly dependent on the Doppler factor 
and its distance from the super-massive black hole in the galac¬ 
tic center, we used the EM cascade argument to determine a min¬ 
imum distance for each Doppler factor value. Finally, for a given 
source size, f var sets an additional lower limit on the Doppler factor 
(d var ). With those additional elements, we arrived at the following 
robust results: (i) the Doppler factor of the emission region should 
be higher than 18-20 for magnetic field strengths < 30 G and 
~day timescale variability; (ii) the y-ray emission region should 
be located outside the BLR, namely at r > IORblr ~ 3 pc; (iii) 
shorter variability timescales, e.g. < 12 hr, push both the minimum 
Doppler factor and distance to even higher values; (iv) the mag¬ 
netic field strength required by the model at pc scale distances is 
stronger than that inferred from observations; and (v) the jet power 
exceeds by at least one order of magnitude the FSRQ accretion 
power. In conclusion, our results disfavour the proton synchrotron 
model for the FSRQ 3C 273. 
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